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Abstract 

We show that forward proton tagging may significantly enlarge the potential of studying the Higgs 
sector at the LHC. We concentrate on Higgs production via central exclusive diffractive processes (CEDP). 
Particular attention is paid to regions in the MSSM parameter space where the partial width of the Higgs 
boson decay into two gluons much exceeds the SM case; here the CEDP are found to have special advantages. 

1 Introduction 



The Higgs sector is (so far) an elusive basic ingredient of the fundamental theory of particle interactions. 
Searches for Higgs bosons, and the study of their properties, are one of the primary goals of the Large Hadron 
Collider (LHC) at CERN, which is scheduled to commence taking data in the year 2007. The conventional 
folklore is that (under reasonable model assumptions) at least one Higgs boson should be discovered at the 
LHC. In particular, if the light Higgs predicted by the Standard Model (SM) exists 1 , it will almost certainly 
be found at the LHC in the first years of running or even maybe before, at the Tevatron. Moreover the LHC 
should provide a complete coverage of the SM Higgs mass range. 

However there is a strong belief that the Standard Model, in its minimal form with a single Higgs, cannot be 
the fundamental theory of particle interactions. Various extended models predict a large diversity of Higgs- like 
bosons with different masses, couplings and even CP-parities. The most elaborated extension of the Standard 
Model is, currently, the Minimal Supersymmetric Standard Model (MSSM) in which there are three neutral (h, 
H and A) and two charged (H + ,H~) Higgs bosons, where h and H are CP-even (m/j < ran) and A is CP-odd. 

lr The theoretical progress in the calculations of the signal and background processes has been enormous in recent years, reaching 
an accuracy of the next-to-next-to-leading QCD corrections for some relevant cross sections, see for example and references 
therein. At the same time, experimentally related issues of the Higgs search physics have been thoroughly investigated, see for 
example 



Just as for the Standard Model, this benchmark SUSY model has been studied in great detail; for a recent review 
see 0. In the MSSM, the properties of the Higgs sector are characterized by the values of two independent 
input parameters, typically chosen to be the pseudoscalar Higgs boson mass, tua, and the ratio, tan/?, of the 
vacuum-expectation-valucs of the two Higgs doublet fields. At tree level, the pseudoscalar A does not couple to 
the gauge bosons and its couplings to down- (up-) type fermions are (inversely) proportional to tan (3. Within 
the MSSM, the mass of the ft,-boson is bounded to to/, ^ 135 GeV (see, for example, 0] and references therein), 
while the experimental 95% CL lower limit for a neutral scalar Higgs is mu — ttia ~ 92 GeV 

Beyond the Standard Model the properties of the neutral Higgs bosons can differ drastically from SM 
expectations. In some extended models, the limit for a neutral Higgs can go down to below 60 GeV. This 
occurs, for instance in the MSSM with explicit CP-violation, where the mass eigenstates of the neutral Higgs 
bosons do not match the CP eigenstates h,H,A, see, for example, |SJ. Further examples are the models with 
extra dimensions where a Higgs can mix with the graviscalar of the Randall-Sundrum scenario {?], which is 
called the radion (see for example [5]). In the latter case, due to the trace anomaly, the gluon-gluon coupling is 
strongly enhanced, which makes this graviscalar especially attractive for searches in gluon-mediated processes, 
assuming that the background issues can be overcome. These extended scenarios would complicate the study 
of the neutral Higgs sector using the conventional (semi)inclusive strategies. 

After the discovery of a Higgs candidate the immediate task will be to establish its quantum numbers, 
to verify the Higgs interpretation of the signal, and to make precision measurements of its properties. The 
separation and identification of different Higgs-like states will be especially challenging. It will be an even more 
delicate goal to establish the nature of a newly-discovered heavy resonance state. For example, how can one 
discriminate between the Higgs of the extended SUSY sector from the graviscalar of the Randall-Sundrum 
scenario (or, even worse, from a mixture of the Higgs and the radion). As was shown in [§], the central exclusive 
diffractive processes (CEDP) at the LHC can play a crucial role in solving these problems, which are so vital 
for the Higgs physics. These processes are of the form 

pp^p + cj) + p, (1) 

where the -I- signs denote the rapidity gaps on either side of the Higgs-like state <fi. They have unique advantages 
as compared to the traditional non-diffractive approaches ^| ^2 ■ These processes allow the properties of tj> to 
be studied in an environment in which there are no secondaries from the underlying events. In particular, if the 
forward protons are tagged, then the mass of the produced central system </> can be measured to high accuracy 
by the missing mass method. Indeed, by observing the forward protons, as well as the <fi — > bb pairs in the 
central detector, one can match two simultaneous measurements of the <fi mass: = Am m i ss i ng and = to^. 
Moreover, proton taggers allow the damaging effects of multiple interactions per bunch crossing (pile-up) to 
be suppressed and hence offer the possibility of studying CEDP processes at higher luminosities [TI]. Thus 
the prospects of the precise mass determination of the Higgs-like states, and even of the direct measurements 
of their widths and 4> — > bb couplings, looks feasible using these processes. A promising option is to use the 
forward proton taggers as a spin-parity analyser JE5J- This may provide valuable additional (and in some cases 
unique) information in the unambiguous identification of a newly discovered state. 

In Section [21 we illustrate the advantages of CEDP, QJ, in exploring the Higgs sector in specific parameter 
ranges of the MSSM. First in Section |2~T1 we consider the so-called "intense-coupling" regime |12U13| where the 
77, WW* and ZZ* Higgs decay modes are suppressed but where, on the other hand, the CEDP cross section 
is enhanced, in comparison with the SM. Next, in Section 2.2, we discuss the decoupling limit (rriA > 2Mz, 
tan/3 > 5) where the light scalar, h, looks very similar to the SM Higgs. In this case the discovery of the heavier 
scalar H is crucial to establish the underlying structure of the electroweak-symmetry-breaking dynamics. If 
mH < 250 GeV, then the H boson may be observed and studied in CEDP. In Section 2.3 we show that CEDP 
may cover the "windows" where, once the h boson is discovered, it is not possible to disentangle the H or A 
bosons by traditional means. In Section 3 we explain how studies of CEDP may help to determine the properties 
of the Higgs-like states, once they have been discovered. 

The plan of the paper is first to present the predictions and then, in Section 4, to explain how the cross 
sections for the CEDP of the <fi CP-even and CP-odd Higgs states are calculated. In this Section we also 
enumerate the uncertainties in the predictions. Simple formulae to approximate the CEDP cross sections for 
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pp — > p + </> + p with cj> — > 66, and the corresponding QCD 66 background, are given in Section 5. These allow 
quick, reliable estimates to be made of the signal-to-background ratio for a Higgs of any mass, and both CP 
parities. 

In Section 6 we focus attention on ways to identify the pseudoscalar boson A. It is unlikely that the pure 
exclusive diffractive production process will have sufficient events, and so possible semi-inclusive signals are 
proposed and studied. 

2 Discovery potential of diffractive processes 

The main aim of this paper is to show how the Higgs discovery potential at the LHC may be enlarged by using 
the unique advantages of forward proton tagging. Special attention is paid to the cases where the partial widths 
for the decay of Higgs bosons into two gluons greatly exceed the SM values. As mentioned above, we shall 
illustrate the advantages of using forward proton tagging to study the Higgs sector by giving predictions for a 
few examples within the MSSM model. 

2.1 The intense coupling regime 

The "intense-coupling regime" ^1 EH is where the masses of all three neutral Higgs bosons are close to each 
other and the value of tan/3 is large. Then the 77, WW* , ZZ* decay modes (which are among the main 
detection modes for the SM Higgs) are strongly suppressed. This is the regime where the variations of all 
MSSM Higgs masses and couplings are very rapid. This region is considered as one of the most challenging for 
the (conventional) Higgs searches at the LHC, see and references therein. On the other hand, here the CEDP 
cross sections are enhanced by more than an order of magnitude (due to the large gg — ► cj> couplings). Therefore 
the expected significance of the CEDP signal becomes quite large. Indeed, this is evident from Fig. 1, which 
shows the cross sections for the CEDP production of h, H, A bosons as functions of their mass for tan/3 = 30 
and 50. The cross sections are computed as described in Section 5, with the widths and properties of the Higgs 
scalar (h, H) and pseudoscalar (A) bosons obtained from the HDECAY code, version 3.0 (TJ|, with all other 
parameters taken from Table 2 of (TJ|; also we take IMODEL = 4, which means the radiative corrections are 
included according to Ref . JH] . 

Let us focus on the main <fi — > 66 decay mode 2 . This mode is well suited for CEDP studies, since a P-even, 
J z = selection rule ^OEIEII su PP resses the QCD 66 background at LO. Indeed, at the LHC, the signal- 
to-background ratio is S/B ~ 3 for a Standard Model Higgs boson of mass m — 120 GeV, if the experimental 
cuts and efficiencies quoted in j 1 1 1 are used. This favourable ratio was obtained since it was argued that proton 
taggers can achieve a missing mass resolution of Am m i ss j ng ~ 1 GeV. Due to the enhancement of the production 
of MSSM scalars, for large tan/3, the QCD background becomes practically negligible. Note that, normally, 
to estimate the statistical significance of the measurement of the signal cross section, we use the conservative 
formula \J S + B for the statistical error. However when we are looking for a new particle, which gives some peak 
on the top of the background, the statistical error may be much lower. It is given just by the fluctuation of the 
background - \[B. On the other hand for a low number of events, it is appropriate to use a Poisson, rather than 
a Gaussian, distribution. Thus the statistical significance to discover a new particle has to be evaluated more 
precisely. For instance, with the cuts and efficiencies quoted in [11], we estimate that the QCD background is B 
= 40(4) events for an integrated luminosity of £ = 300 fb" 1 (30 fb" 1 ), if we take M = 120 GeV and AM = 1 
GeV. Thus, to obtain a confidence coefficient of 5.7 x 10~ 7 (equivalent to a statistical significance of 5a for a 
Gaussian distribution |2]), it is sufficient for the cross section of a Higgs signal, times the 66 branching fraction, 
to satisfy 

Br(66) • cr > 0.7 fb (2.7 fb) (2) 

2 The exclusive diffractive studies in Ref. 1111 addressed mainly the bb decay mode. The use of the tt decay mode requires 
an evaluation of the pp — » p + tt + p background, especially of the possibility of misidentifying gluon jets as t's in the CEDP 
environment. Currently, for inclusive processes, the probability, P(g/r), to misidentify a gluon as a t is estimated to be about 
0.002 [T7| 



3 



Central exclusive diffractive production 
o Br(h/H->bb) (fb) a Br(h/H^>bb) (fb) 



h:H 



tanp = 30 



SM 



h:H 



tanP = 50 



SM 



10 u 1 ^ 1 1 1 1 1 10 

100 150 200 250 300 100 150 200 250 300 



m h/H (GeV) 



m h/H (GeV) 



a Br (fb) 



tanp = 30 




cBr (fb) 




m A (GeV) 



150 175 200 

m A (GeV) 



Figure 1: The cross sections, times the appropriate bb and t + t~ branching fractions, predicted for central 
exclusive diffractive production of h(0 + ), H(0 + ) and ^4(0~) MSSM Higgs bosons (for tan/3 = 30 and 50) at the 
LHC, obtained using the MRST99 ^3] gluon distribution. The dotted curve in the upper plots shows the cross 
section for the production of a SM Higgs boson. The vertical line separates the mass regime of light h(0 + ) and 
heavy if(0 + ) bosons. The curves were computed using version 3.0 of the HDECAY code JS], with all other 
parameters taken from Table 2 of [TJ]; the radiative corrections were included according to Ref. |16|. 



if we take the integrated luminosity to be C — 300 fb 1 (30 fb 1 ). 

As can be seen from Fig. 1, for MSSM with tan/? = 50 the expected CEDP cross section times branching 
fraction, Br (if — > bbjan, is greater than 0.7fb for masses up to irijj ~ 250 GeV. The situation is worse for 
pseudoscalar, A, production, since the cross section is suppressed by the P-even selection rule. Thus the CEDP 
filters out pseudoscalar production, which allows the possibility to study pure if production, see Fig. 1. This 
may be also useful in the decoupling limit, to which we now turn. 

2.2 The decoupling limit 

For rriA > 2m z and tan/3 > 5, there is a wide domain in MSSM parameter space where the light scalar 
h becomes indistinguishable from the SM Higgs, and the other two neutral Higgs states are approximately 
degenerate in mass. This is the so-called decoupling limit 22 . Moreover in many theories with a non-minimal 
Higgs sector, there are significant regions of parameter space that approximate the decoupling limit; see, for 
instance, [3] . In this regime the discovery of the heavier non-minimal Higgs scalar if is crucial for establishing 
the underlying structure of the electroweak-symmetry-breaking dynamics. Here, forward proton tagging can 
play an important role in searching for (at least) the if -boson, if it is not too heavy (m# < 250 GeV). For large 
values of tan/3 the decoupling regime essentially starts at ~ 170 GeV. As can be seen in Fig. 1, the cross 
section is still sufficiently large to ensure the observation of the if boson up to tcia — 250 GeV. To illustrate 
the moderate tan/3 region we consider the case tua = 185 GeV and tan/3 = 7. Then the CEDP cross section is 

Br(if ->• bb)a H = 0.17 fb. (3) 
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Thus for an integrated luminosity of C = 300 fb _ we have about 8.5 events, after applying the experimental 
cuts and efficiencies discussed in Ref. jllj . The corresponding background is about 2.6 events. Using again a 
Poisson distribution, as we did in the case of ©, then the expected statistical significance to observe such an 
H boson is about 4ct. So detection may be feasible. 

The possibility to use exclusive diffractive processes to explore larger masses will depend on various experiment- 
related factors. In particular, on the prospects to achieve better mass resolution, Am, at higher mass, m. 

2.3 The 'window' or 'hole' regions 

There are regions of MSSM parameter space where, after the h boson is discovered, it will be difficult to identify 
the heavy H or A bosons by conventional non-diffractive processes. Can the CEDP processes help here? First 
consider the 'notorious hole' in parameter space (sometimes called the "LHC wedge"), defined by tua ^ 200 GeV 
and tan/? ~ 4-7 [3]. In this region, only the lightest Higgs h can be discovered even with integrated luminosity 
of the order of 300 fb" 1 . Its properties are nearly indistinguishable from those of a SM Higgs. Unfortunately 
here, even for H, the predicted CEDP cross section is too low, assuming the efficiencies and cuts of Ref. 
see the heavy continuous curve in Fig. 2. 
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Figure 2: The cross sections for the central exclusive diffractive production of h and H bosons, pp — > p+(h, H)+p, 
times the appropriate bb branching fraction, as a function of their mass, in two regions of MSSM parameter 
space where observation of the H using non-diffractive processes is difficult. Note that in the first region the 
curve changes from dotted to continuous when the mass of the ^4(0~) bosons reaches 200 GeV; and that the 
value of nih is practically unchanged for values of tua in the region rriA > 200 GeV. 



However, with the present understanding of the capabilities of the detectors it is plausible that similar 
windows will occur in the interval 130 < tua ^170 GeV for tan (3 ~4-6. 3 That is, it may be impossible to identify 
the H scalar by conventional processes at the ha confidence level with 300 fb -1 of combined ATLAS+CMS 
luminosity. From the predictions shown by the dashed line in Fig. 2, we see that the part of this 'hole', up to 

3 We are grateful to Sasha Nikitenko for pointing out the importance of the additional (nonconventional) studies in this region. 
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inn ~ 160 GcV, may be covered by CEDP. For illustration, consider the case itih = 150 GeV and tan/3 = 4.5. 
With an integrated luminosity of C — 300 fb -1 , the CEDP cross section of 

Br(i? -> 66)cr ff = 0.8 fb (4) 

gives about 40 events, after applying the experimental cuts and efficiencies discussed in Ref. The corre- 

sponding background is about 10 events. Thus, if we use a Poisson distribution, as we did in the case of J2J, 
then the expected statistical significance to observe such an H boson will exceed 10er. 

2.4 Exotic scenarios 

Recently some interesting applications of the CEDP to studies of the more exotic Higgs scenarios, such as 
models with CP violation and searches for the radions, have been discussed |23|. However the feasibility of 
applying these studies will require further critical discussion of the background issues. 

3 Potential of diffractive processes after a Higgs discovery 

As was discussed in Refs. ^JE], the use of CEDP, with forward proton tagging, can provide a competitive way 
to detect a Higgs boson. This approach appears to have special advantages for searches of the CP-even Higgses 
of MSSM in the intense-coupling regime. Indeed, the CEDP cross section can exceed the SM result by more 
than an order of magnitude, see Fig. 1. Once a Higgs- like signal is observed, it will be extremely important 
to verify that it is indeed a Higgs boson. Here we describe where forward proton tagging can provide valuable 
help to establish the nature of the newly discovered state and to verify whether it can be interpreted as a Higgs 
boson. We also emphasize that CEDP are well suited to study specific cases which pose difficulties for probing 
the Higgs sector at the LHC. 

It is generally believed that a detailed coverage of the Higgs sector will require complementary research 
programmes at the LHC and a Linear e + e~ collider. Interestingly, in some sense, the implementation of 
forward proton taggers allows some studies which otherwise would have to await the construction of a Linear 
e + e~ collider. Moreover, in addition to the other advantages, forward proton tagging offers the possibility to 
study New Physics in the photon-photon collisions, see for example |24l I25ll2l)lll0| . 

As mentioned above, if a candidate Higgs signal is detected it will be a challenging task to prove its Higgs 
identity. Unlike the conventional inclusive approaches, the very fact of seeing the resonance state in CEDP 
automatically implies that the newly discovered state has the following fundamental properties. It must have zero 
electric charge and be a colour singlet. Furthermore, assuming P and C conservation, the dominantly produced 
state has a positive natural parity, P = (— 1) J and even CP. Recall that the installation of forward proton 
taggers also allows the attractive possibility of a spin-parity analysis. This may provide valuable additional 
(and in some cases unique) leverage in establishing the origin of the discovered candidate state. In particular, 
assuming CP conservation, the CEDP allow the _ , 1 _ , 1 + states to be filtered out, leaving only an ambiguity 
between the ++ and 2 ++ states. Though without further efforts the 2 ++ state cannot be ruled out, this would 
not be the most likely choice. 

As discussed in [5], studying of the azimuthal correlations of the outgoing protons can allow further spin- 
parity analysis. In particular, it may be possible to isolate the 0~ state. The azimuthal distribution distinguishes 
between the production of scalar and pseudoscalar particles [§]• Note that with the forward protons we can 
determine the CP-properties of the Higgs boson irrespective of the decay mode. Moreover, CEDP allow the 
observation of the interference effects between the CP-even and CP-odd gg — > cj) transitions (see |27|L 

Following the discovery, and the determination of the basic quantum numbers, of the new state, the next 
step is to check whether the coupling properties match those expected for the Higgs boson. For this, we need 
to detect the signal in at least two different channels. The CEDP mechanism provides an excellent opportunity 

4 to detect the bb mode. The other decay channel which is appropriate for CEDP studies is the rr decay mode 
(or WW*, or /u/i, if the event rate is sufficient). These measurements allow a non-trivial check of the theoretical 

4 Recall that for more inclusive processes the signal is overwhelmed by the QCD bb background. 
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models. For example, the precise measurement of the partial decay widths into tau and bottom quark pairs can 
provide a way to discriminate between supersymmetric and non-supersymmetric Higgs sectors, and to probe 
CP-violating effects (for recent papers see, for instance, and references therein). 

For illustration, consider the following topical example, namely the intense coupling limit of MSSM. As 
mentioned before, this regime is especially difficult for conventional inclusive studies of the MSSM Higgs |12U13| . 
On the other hand, with forward proton taggers, CEDP offer an ideal probe. 

For numerical purposes, we choose the same parameters as in |13| . in particular, tan/3 = 30, together with 
results for tan (3 — 50. The salient features of this regime are that there is almost a mass degeneracy of the neutral 
Higgs states, and their total widths can be quite large (due to the tan 2 /3 increase) and reach up to 1-2 GeV. 
This can be seen from Fig. 3, where we plot the mass differences AM = ttia — and AM = mn — m^, and 
the widths of the bosons, as a function of the pseudoscalar mass tua- It is especially difficult to disentangle the 
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Figure 3: The upper plots show the differences AM — tua — rrih (continuous curve) and AM = mn — ttla 
(dashed curve) as a function of tua for tan (3 = 30 and 50 respectively. The lower plots show the total widths 
of the MSSM neutral high bosons. 



Higgs bosons in the region around tua ~ 130 GeV, where all three neutral Higgs states practically overlap |13j . 
Since the traditional non-diffractive approaches do not, with the exception of the 77 and fifi modes, provide a 
mass resolution better than 10-20 GeV, all three Higgs bosons will appear as one resonance 5 . 

As an example, we give the expected number of CEDP signal and background events for C — 30 fb _1 for 
rriA — 130 GeV, using the efficiencies and angular cuts adopted in Ref. (TJ. For tan/3 = 30(50), the number of 
the A ^ bb signal events, Sa — 0.35(1), while the QCD background is B = 2.5 for AM = 1 GeV. For this choice 
of parameters, rrih — 122.7(124.4) GeV, and the number of the h — > bb events is Sh = 30(71) with background 
B = 3.5(3.2). The corrresponding numbers for the H boson are m H = 134.2(133.5) GeV, S H = 47(124) and 

5 Recall that in this regime the 77 decay mode is hopeless. Moreover, the vector-boson-fusion Higgs production, with the bb or 
tt decays, does not allow the separation of the h and H peaks due to the poor invariant mass resolution. As shown in 1131 . only the 
detection of the rare dimuon Higgs decay mode may provide a chance to separate the Higgs states in inclusive production. Even 
this would require that the Higgs mass splitting exceeds at 3-5 GeV. 
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B = 2.1. 

An immediate advantage of CEDP, for studying this troublesome region, is that the A contribution is 
strongly suppressed, while the h and H states can be well separated (m# — m,h — 10 GeV) given the anticipated 
experimental mass resolution of AM ~ 1 GeV JT], see Fig. 4. Note that precision measurements of the two 




m A (GeV) m A (GeV) 



Figure 4: The mass bands ± T for neutral MSSM Higgs bosons as a function of rriA, for different values of 
tan p. A vertical slice through these plots, at a given value of tua, indicates the position of the MSSM resonance 
Higgs peaks. The upper right hand plot shows that the h and H bosons are clearly identifiable for tan (3 = 30, if 
A(0~) production is suppressed. The lower plots show how the sensitivity of the widths, to variations of tan/3, 
will change the profile of the peaks. This way of presentation was motivated by Ref. [131 • Note that the widths 
of the bosons for tan/3 less than about 15 become barely visible on this type of plot. 



CP-even Higgs states may allow a direct determination of the basic MSSM parameter, tan/3. 

Moreover, the forward tagging approach can provide a direct measurement of the width 6 of the h (for 
rnu ^5 120 GeV) and the iJ-boson (for run ^ 130 GeV). Outside the narrow range m,A — 130 ± 5 GeV, the 
widths of the h and H are quite different (one is always much narrower than the other), see Fig. 3. Recall 
that the width measurement can play an important role in the understanding of Higgs dynamics. It would be 
instructive to observe this phenomenon experimentally. 

As discussed in for tan/3 = 30 the central exclusive signal should be still accessible at the LHC up to an 
H mass about 250 GeV. For instance, for m H = 210 GeV, and LHC luminosity 30 ftT 1 (300 fb" 1 ), about 20 
(200) H — > bb events are produced. If the experimental cuts and efficiencies quoted in DKMOR are imposed, 
then the signal is depleted by about a factor of 6. This leaves 3 (30) observable events, with background of 
about 0.1 (1) events. 

6 Note that, strictly speaking, when the Higgs width, T, becomes comparable to the mass resolution, Am, the zero-width 
approximation is not valid any more. Therefore in calculations of the cross sections the mass resolution and large-width effects 
should be taken into account. 
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(a) (b) 




Figure 5: (a) The QCD diagram for double-diffractive exclusive production of a Higgs boson h, pp — > p + h + p, 
where the gluons of the hard subprocess gg — > h are colour screened by the second i-channel gluon. (b) The 
rescattering or absorptive corrections to pp — > p + h + p, where the shaded region represents the soft pp 
rescattering corrections, leading to the suppression factor S 2 . 



4 Calculation of the ± Higgs cross sections and their uncertainties 

To calculate the cross section for the central exclusive diffractive production of ^(O*) Higgs bosons 7 we use the 
formalism of Refs. [5511201 EH- The amplitudes are described by the diagram shown in Fig. [3a), where the hard 
subprocesses gg — > ft-(0 ± ) are initiated by gluon-gluon fusion and where the second i-channel gluon is needed 
to screen the colour flow across the rapidity gap intervals. Ignoring, for the moment, the screening corrections 
of Fig. ISJb), the Born amplitudes are of the form [§] 

T h = An 2 [ d y^ V ^ ^-zfg(xi,x' 1 ,Q 2 1 ,fi 2 - 1 t 1 )f g (x2,x' 2 ,Ql,fi 2 - 1 t 2 ), (5) 

J QAQ± - pi±) 2 {Q± + P2±) 2 

where the gg — > h^^) subprocesses are specified by [291 I30| 



A% = (6) 
with the NLO K factor K ~ 1.5 |3THTD] . and by the factors 



V h (o+) = (Qjl -fix) ■ (Q_L +P2±) 

Vh(o~) = (iQ± ~ Pi-l) x (Q± + p2±fj ■ no ■ 



(7) 



Here Hq is a unit vector in the beam direction. The / g 's are the skewed unintegrated gluon densities of the 
proton at the hard scale ^i, taken to be with 



Qi = min|<3x, \{Qx — Pu 
Q 2 = min Iq±, \(Q± + p 2j 

The longitudinal momentum fractions carried by the gluons satisfy 



(8) 



Qx ^«(x~^)«l (9) 



where, for the LHC, with y/s = 14 TeV, we have x ~ 0.01, while x' ~ 10 -4 . Below, we assume factorization of 
the unintegrated distributions, 

f g (x, x', Q 2 , fj, 2 : t) = f g (x, x', Q 2 ^ 2 )F N (t), (10) 



7 For convenience of presentation in this section we denote the pseudoscalar boson by h(0 ) rather than A(0 ). 
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where we parameterize the form factor of the proton vertex by the form Fjv(i) = exp(6i) with 6=2 GeV" 2 . In 
the domain specified by the skewed unintegrated densities are given in terms of the conventional (integrated) 
densities g(x,Q 2 ). To single log accuracy, we have |32] 8 . 



f g (x,x' ',Ql,n 2 ) 



Rc 



d 



(^/T{Qi,n) xg{x,Ql 



(11) 



where T is the usual Sudakov form factor which ensures that the gluon remains untouched in the evolution up 
to the hard scale fx, so that the rapidity gaps survive. This Sudakov factor T is the result of resumming the 
virtual contributions in the DGLAP evolution. It is given by 



,~ s ( f' a s(k 2 ) dkf f 1 '* 



1 



dz 



(12) 



Here we wish to go beyond the collinear approximation and in the T factor to resum, not just the single collinear 
logarithms, but the single soft log 1/(1 — z) terms as well. Thus we consider the region of large-angle soft gluon 
emission, and explicitly calculate the one-loop vertex diagram. In particular, we do not approximate sin 9 by 6. 
Then we adjust the upper limit of the z integration in l|12|) to reproduce the complete one- loop result 9 . We find 



k t + 0.62m h 



(13) 



The square root in (|llf> arises because the bremsstrahlung survival probability T is only relevant to hard gluons. 
R g is the ratio of the skewed x' <C x integrated distribution to the conventional diagonal density g(x,Q 2 ). For 
x -C 1 it is completely determined |33| , with an uncertainty in R g of the order of 4x 2 . The apparent infrared 
divergence of JHJ) is nullified 10 for h(0 + ) production by the Sudakov factors embodied in the gluon densities f g . 
However, as discussed in Ref. [H], the amplitude for h(0~) production is much more sensitive to the infrared 
contribution. Indeed let us consider the case of small pi± of the outgoing protons. Then, from J2J, we see 

Pi±P2± (since the linear contribution in Q± vanishes after the angular 



that V y 



h(o+) 



f ± , whereas V h{0 -) 
i2n . /r>4 



integration). Thus the d 2 Q±/Q 4 L integration for h(Q + ) is replaced by pi±P2±d 2 Q±/Q e L for h(0 ), and now the 
Sudakov suppression is not enough to prevent a significant contribution from the Q\ < 1 GeV 2 domain. 



4.1 Uncertainties in predictions of the cross sections 

To estimate the uncertainty in the predictions for the h (0) exclusive diffractive cross sections we first quantify 
the above uncertainty arising from the infrared region, where the gluon distribution is not well known. As an 
example, consider h(0~) and h(0 + ) production at the LHC, for rrih — 120 GeV and /i = m^/2, using different 
treatments of the infrared region. We perform calculations using MRST99 ^3] and CTEQ6M partons 
respectively with the very low Q gluon frozen at its value at Qi,2 — 1-3 GeV. Then we integrate down in Q± 
until Qi.2 are close to Aqcd, where the contribution vanishes due to the presence of the T-factor. This will 
slightly overestimate the cross sections, as for x ~ 0.01 the gluon density has a positive anomalous dimension 
(xg ~ (Q 2 ) 7 with 7 > 0) and decreases with decreasing Q 2 . A lower extreme is to remove the contribution 
below Qi.2 = 1.3 GeV entirely. Even with this extreme choice, the + cross section is not changed greatly; it 
is depleted by about 20%. On the other hand, as anticipated, for 0~ production, the infrared region is much 
more important and the cut reduces the cross section by a factor of 5. 

Another uncertainty is the choice of factorization scale p.. Note that in comparison with previous calcu- 
lations jl()| . which were done in the limit of proton transverse momenta, pi t 2± <C Q±, now we include the 
explicit p^-dependence in the Q±-\oop integral of JSJ. We resum the 'soft' gluon logarithms, In 1/(1 — z), in 

8 In the actual computations we use a more precise form as given by eq.(26) of Ref. 1321 

9 An explicit computation of the one-loop vertex contribution from the region of gluon transverse momentum <C m h/2 
gives (N c us /it) [0.212 + ln(m^/2fc^) — ll/12]dfc^/fc^. The last number in the square brackets corresponds to the standard non- 
logarithmic contribution related to the zP gg (z) term in 1121 . while the first number stems from the exact treatment of large-angle 
gluon radiation. This result can be incorporated in 11 21 by simply choosing A as in 1131 . 

10 In addition, at LHC energies, the effective anomalous dimension of the gluon gives an extra suppression of the contribution 
from the low Q± domain 30 . 
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the T-factor. So now the T-factor includes both the soft and collinear single logarithms. The only uncertainty 
is the non-logarithmic NLO contribution. This may be modelled by changing the factorization scale, //, which 
fixes the maximal k t of the gluon in the NLO loop correction. As the default we have used /z = m^/2; that is 
the largest k t allowed in the process with total energy m^. Choosing a lower scale fi = rrih/A would enlarge the 
cross sections by about 30%. Recall that the general calculation is performed in the collinear approximation. 
However, when we take account of large-angle soft gluon emission in the form factor, that is in the calculation of 
the virtual- loop correction, as in (|12|) and i|13|) . we go beyond the collinear approximation. The corresponding 
real contribution in is represented by R g and the derivative of the \/T factor. Strictly speaking, for the 
case of large-angle emission there will be some small deviation from this simplified formula. This gives an extra 
uncertainty of the order of 10%. 

Next there is some uncertainty in the gluon distribution itself. To evaluate this, we compare predictions 
obtained using CTEQ6M gU, MRST99 [H] and MRST02 [35] partons. For 0+ production at the LHC, with 
rrih = 120 GeV and \x = m^/2, we find that the effective gluon-gluon luminosity, before screening, is 



dC 



dyd In m 2 



= (2.2, 1.7, 1.45) x 10~ 2 (14) 



respectively. This spread of values arises because the CTEQ gluon is 7% higher, and the MRST02 gluon 4% 
lower, than the default MRST99 gluon, in the relevant kinematic region. The sensitivity to the gluon arises 
because the central exclusive diffractive cross section is proportional to the 4th power of the gluon. For 0~ 
production, the corresponding numbers are (4.2, 2.7, 1.7) x 10~ 5 . Up to now, we have discussed the effective 
gluon-gluon luminosity. However, NNLO corrections may occur in the gg —> h fusion vertex. These give an 
extra uncertainty of ±20%. Note that we have already accounted for the NLO corrections for this vertex |10|. 

Furthermore, we need to consider the soft rescattering which leads to a rather small probability, S 2 = 0.026, 
that the rapidity gaps survive the soft pp interaction, see Fig.OJb). From the analysis of all soft pp data 
we estimate the accuracy of the prediction for S 2 is ±50%. One check of the eikonal model calculations of S 2 
is the estimate of the diffractive dijet production rate measured by the CDF collaboration [37] at the Tevatron. 
The rate, when calculated using factorization and the diffractive structure functions obtained from HERA data, 
lies about a factor of 10 above the CDF data. However, when rescattering corrections are included, and the 
survival probabilities computed, remarkably good agreement with the CDF measurements is obtained |38| . 

Combining together all these sources of error we find that the prediction for the + cross section is uncertain 
to a factor of almost 2.5, that is up to almost 2.5, and down almost to 1/2.5, times the default value. 11 On 
the other hand, 0~ production is uncertain by this factor just from the first (infrared) source of error, with the 
remaining errors contributing almost another factor of 2.5. 



5 Simplified formulae for the CEDP signal and background 

Here we give simple approximate formulae for the cross sections for the central exclusive diffractive production 
of /i(0 ± ) bosons, pp — > p + h + p, and for the QCD bb background as a function of mass. As was shown in 
Ref. ^Ojj the cross section may be written as the product of the effective gluon-gluon luminosity C and the 
cross section for the hard subprocess a 

CCEDP = C&. (15) 

At the LHC energy, the luminosity, integrated over rapidity, may be approximated, to an accuracy of better 
than 10% in the interval 25 < m < 250 GeV, by 



dC 



dim 



dC 



dyd In m? 



dy = < 



9.1 x 10 5 (S 2 ) 0+ 
(16±to) 3 - 3 

1.94 x 10 3 (S 2 ) - 
(0.5 ±m) 3 
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for h(0 + ) production 



for h(0 ) production 



(16) 



11 For example, we predict the cross section for the exclusive diffractive production of a Standard Model Higgs at the LHC, with 
rrih = 120 GeV, to be 2.2 fb with an uncertainty given by the range 0.9—5.5 fb. 
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where the mass m is in GeV. The factors (S 2 ), which are the probability that the rapidity gaps survive the soft 
rescattering, are 0.026 and 0.087 for h(0 + ) and h(0~) respectively. The 0~ amplitude contains the kinematic 
factor pi± x P2j_, which implies that in impact parameter space the Fourier transform contains the factor 
b%± x &2± EH- Thus the 0~ amplitude tends to populate larger impact parameters, where the suppression caused 
by soft rescattering is less effective, than does the + amplitude. The gg — > h subprocesses cross section is, in 
the zero width approximation, 

. 2n 2 T(h -> gg) ( m 2 

a h = = oil j 

m A h V m h 

where Gf is the Fermi constant. 

The bb background processes were discussed in detail in Ref. m for a Standard Model h — ► bb signal, with 
rrih = 120 GeV. There were four appreciable sources of bb background, each with a background-to-signal ratio 
in the range B/S = 0.06-0.08, after the appropriate cuts. In brief, these give contributions to the background 
subprocess cross sections of the following approximate forms 

(i) gg pp — * gg (bb mimicked by gg): a ~ a|/m 2 , 

(ii) gg pp — ► bb (J z — 2 admixture): a ~ a 2 /m 2 , 

(hi) gg pp — * bb (rn& 7^ 0, J z = contribution): a ~ a 2 s m 2 /m i , 
(iv) NLO gg pp — > 66g contribution: <r ~ a|/m 2 , 

where the renormalisation scale is specified by as(m 2 /4). The PP superscript is to denote that each incoming 
gluon to the hard subprocess belongs to a colour-singlet t-channel state, see Fig. [S] Note that all background 
contributions go like a ~ l/m 2 , except for (in). In addition, for fixed missing mass resolution Am, the 
background contains a factor Am/m. Thus for large m > 100 GeV the background cross section behaves 
roughly as 

a B = Ca B - Am/m 6 , (18) 

while the Higgs signal behaves as 

o h = Co h - l/m 3 . (19) 

To obtain a more precise behaviour one can use (|16fl . 

In Ref. we found a signal-to-background ratio S/B ~ 3 for a m^ = 120 GeV Standard Model Higgs 
boson. The ratio may, however, be appreciably enhanced for SUSY Higgs with large m, and/or large tan/3. It 
is interesting to note that, in the region where the background is larger than the signal (B > S), the statistical 
significance of the signal S/^JS + B ~ m~ 3 / y/ Am/m 6 is practically independent of mass, assuming that the 
experimental resolution Am is mass independent. Moreover, for a low mass, where bb pairs in the J z — state, 
(iii), give the dominant contribution to the background, the statistical significance increases with mass, as here 
B - Am/m 8 . 

Finally, it is important to note that there are still various open questions which must be addressed by 
experimentalists in order to establish the feasibility of CEDP for Higgs searches at the LHC. One of the most 
challenging issues concerns triggering. Assuming that the trigger is based on the central detector only, the main 
task is to reduce the rate of triggering jets in order not to reach saturation of the trigger. We present below 
some numbers which may be useful when considering this issue. 

Based on the results of Ref. ^UJi we see that the cross section for the CEDP of dijets with Et > 30 GeV is 
200 pb. If we assume that the proton tagging efficiency is 0.6, then this yields 0.12 events/sec for a luminosity 
of 10 33 cm _2 sec _1 . At the same time, the inclusive cross section for production of a pair of jets, each with 
Et > 30 GeV is about 50 fib. This leads to 50,000 events/sec, which greatly exceeds the trigger saturation 
limit. For CMS, for example, this latter limit is 100 events/sec [35]. Therefore it appears necessary to devise 
additional topological requirements in the level- 1 trigger for the exclusive signal, such as making use of the 
presence of the rapidity gaps. 



a: 



(17) 
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6 Identifying the pseudoscalar boson, A 



If the exclusive diffractive cross sections for scalar and pseudoscalar Higgs production were comparable, it would 
be possible to separate them readily by the missing mass scan, and by the study of the azimuthal correlations 
between the outgoing proton momenta. However, we have seen that the cross section for pseudoscalar Higgs 
exclusive production is, unfortunately, strongly suppressed by the P-even selection rule, which holds for CEDP 
at LO. For values of tan/3 ~ 10-15, the separation between the H/h and A bosons is much larger than their 
widths, see Fig. 4(d). Hence it might be just possible to observe the pseudoscalar in CEDP. For example, for 
tan/3 = 15 and wia — 95 GeV, the mass separation, 3.6 GeV, between h and A is about 8 times larger than the 
width T/j. The cross section 

Br(j4 -► 66) • a A ~ 0.15 fb, (20) 

when allowing for the large uncertainties in a a, could be just sufficient to bring the process to the edge of 
observability. 

On the other hand, one may consider the semi-exclusive process 

pp^p + j<l>j+p, (21) 

where the Higgs bosons, (f> — H, A, are accompanied by two gluon jets with transverse momenta p± satisfying 
m 4, P± ^ 1 GeV. The emission of soft gluons, with energies (in the </> rest frame) E <C does not affect 

the P-even selection rule. Thus in this region the scalar H boson will be dominantly produced in process (|2ip. 
On the other hand, in the domain 

x t ee 2E i /m 4> ~ 1 (22) 

we lose the selection rule and have comparable H and A production. Moreover the azimuthal angular distri- 
bution 12 for the gluon jets is different for the H and A cases. In particular, for gluons with X\ = X2 = 1 
accompanying A production, the azimuthal distribution is practically flat (at LO), whereas for H production 
it has the form 3 + 2cos(2</>). This observation may be used to distinguish between H and A production, and 
to study better the properties of the A boson. 

Of course, we have to pay a price, of about a factor 1/20, to produce two extra jets with p± in the range 
3> p± 3> 1 GeV, in a rapidity interval with Alnxj ~ 1. However, this is less than the suppression of A 
production by the P-even selection rule. For tan (3 > 30 and rriA < 200 GeV the cross section still exceeds 1 fb. 

To gain more insight, it is informative to return to the original central exclusive diffractive process, which 
we sketch in Fig.|Sfa). In the equivalent gluon approximation, the polarization vectors q of the gluons, carrying 
momentum fractions Xi, are aligned along the corresponding gluon transverse momenta (Q— Pi)x an d (Q+P2)±- 
For very forward protons with pi^± <C this leads to a correlation between the polarization vectors, namely 
e*i is parallel to e*2- On the other hand, the <f> production vertices behave as ^(0 + ) ~ ei • e*2 and V(0~) ~ e\ X e*2- 
Therefore the exclusive diffractive production of the pseudoscalar A boson is strongly suppressed 13 . 

To change the i-channel gluon polarization we need particle emission. The best possibility would be to emit 
either a rf or a r\ c 0~ meson. Unfortunately the corresponding cross sections are too low. We lose more than 
a factor of 10~ 4 , which is a stronger suppression than that coming from the P-even, J z = selection rule. 
For this reason we considered the possibility of emitting a (colour singlet) pair of gluons, that is process (|21|l . 
For 'soft' (xi « 1) emission, the dominant contribution comes from the term in the triple-gluon vertex which 
contains g vv , and so does not change the gluon polarization vector. Here the indices v and v' correspond to the 
t-channel gluon polarization before and after the emission of the s-channel gluon. However for large Xi > 1 the 
major contribution comes from the term with g v,i , where fi is the s-channel gluon polarization index. In this 
case after the emission, the new i-channel polarization vector is directed along the new transverse momentum 

12 The possibility to discriminate the CP-parity of different Higgs states by measuring the azimuthal correlations of the two tagged 
gluon jets in inclusive <j>gg production was mentioned in Ref. 1401 . However, the proof of the feasibility of such an approach in 
non-diffractive (frgg production requires further study of the possible dilution of the effect caused by parton showers (see also 1411 ). 
An exclusive environment of the jets may be needed. 

13 However, when the polarisation vectors are parallel to Q±, it follows, after integration over the azimuthal angle of Qj_, that 
the LO QCD production of bb pair is also suppressed by a factor (m b / Ex) 2 ; which results in a very low QCD bb background to the 
CEDP. 
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Q' ± ~ —p±. Thus for Xi ~ 1 the cross sections for iJ and A production are comparable, but now the LO QCD 
background is no longer suppressed. 

Maybe the best chance to identify the ^4(0~) boson is to observe the double-diffractive inclusive process 

pp -» X + <f> + Y, (23) 

where both protons are destroyed. In terms of i-channel gluon polarization, this is equivalent to the limit x, 3> 1 
in the previous reaction (|21fl . Process (|23[1 has the advantage of a much larger cross section 14 . However, just as 
for the above semi-exclusive process, we do not have the J z = selection rule to suppress the bb background, 
nor do we have the possibility of the good missing mass resolution 15 . On the other hand, the gg pp luminosity 
is more than order of magnitude larger than for the pure exclusive case. For example, for double-diffractive 
inclusive production, with the rapidity gaps Ai] > 3, the luminosity is 20 times larger than that for the exclusive 
diffractive production of a Higgs boson with mass mu = 120 GeV (see for details J0|)- So, even for the rr 
decay mode we expect a cross section, <r- mc i, of about 20 fb for A production in the MSSM with tan/3 = 30 and 
rriA — 120 GeV. This looks promising, provided that the probability for gluon misidentification as a r is less 
than 1/150, which appears to be feasible |17U42| . 

Despite the QCD background, the observation of the A boson also appears to be feasible for the bb decay 
mode, since then a ~ 200 fb which gives a signal-to-background ratio of S/B = 1/24. We have evaluated the 
QCD bb background using the LO formula for the gg pp — ► bb cross section (without any K factor) with the 
polar angle cut 60° < 8 < 120° as was done in |11|. and assumed a mass resolution of AM = 10 GeV. The 
efficiency of tagging the b jets was taken to be about 0.6 This results in a statistical significance 16 of 

the signal for the A boson of S/y/S + B ~ 8.5er, for an integrated luminosity of 30 fb -1 . Moreover one can 
measure the azimuthal angle ip between the transverse energy Et flows in the forward and backward regions, 
that is between the total transverse momenta of the X and Y systems. The angular dependence is driven by 
the momenta transferred across the rapidity gaps, and has the form 

cos 2 <^ for#(0+) 
and sin 2 ip for A(0~), 

whereas for the QCD bb background the ip dependence is flat (if we neglect the quark mass (mj, <§; M^/2) and 
average over the 6-jet momentum direction; see Eq.(15) of 0Hj for Q = 0). Thus we may suppress the H(0 + ) 
signal by selecting events with <p close to ir/2. We emphasize that here we deal with rather large Et jets, 
Et ^5 — 20 GeV, and so the angular dependences are largely insensitive to the soft rescatterings. Furthermore, 
instead of measuring the pr of the jets (as in Ref. 40 ), for process (|23|l it is enough to measure the total 
transverse energy Et in the forward (and backward) regions of proton dissociation. 

Recall that the effective gg pp luminosity for the inclusive process H23JI was evaluated using the leading 
order formula. Thus we cannot exclude rather large NLO corrections, particularly for the non-forward BFKL 
amplitude. On the other hand, here, we do not enter the infrared domain. Moreover, with "skewed" kinematics, 
the NLO BFKL corrections are expected to be much smaller than in forward case. So we may expect that the 
uncertainty of the prediction is about factor of 3 - 4 (or even better) . 



7 Conclusion 

The Central Exclusive Diffractive Processes promise a rich physics menu for studying the mechanism of elec- 
troweak symmetry breaking, and the detailed properties of the Higgs sector. We have seen, within MSSM, that 
the expected CEDP cross sections are large enough, especially for large tan (3, both in the intense coupling and 

14 Note that instead of the factors 1/fe, which came from the tj = — p?. integrals limited by the proton form factors, as in 1101 
(see also eq.(8) in Ref. 1101 ). here we have logarithmic dti/ti integrals over the intervals <C |t;| <C m?. Moreover, the Sudakov 
supression (which arises from the requirement that no extra relatively 'hard' gluons are emitted with transverse momenta kj_ in 
the interval \Q± ipi±\ «C kj_ <SC m^/2) is much weaker than for the exclusive process due to the larger values of Pix, see 1101 1291 . 
Therefore the effective gg FP luminosity for double-diffractive inclusive production is much larger than in the CEDP case. 

15 Also, in this case, pile-up can cause problems. 

16 If we include a K factor, K ~ 1.5, for the background then the statistical significance is about 7a. 
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the decoupling regimes. Thus CEDP offers a way to cover those regions of MSSM parameter space which may 
be hard to access with the conventional (semi) inclusive approaches. This considerably extends the physics 
potential of the LHC and may provide studies which are complementary, both to the traditional non-diffractive 
approaches at the LHC, and to the physics program of a future Linear e + e~ collider. To some extent this 
approach gives information which otherwise may have to await a Linear e + e~ collider. 
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